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ABSTRACT 


Aims. We investigate the nature of dissipative instability at the boundary (seen here as tangential discontinuity) between 
the viscous corona and the partially ionised prominence plasma in the incompressible limit. The importance of the partial 
ionisation is investigated in terms of the ionisation fraction. 

Methods. Matching the solutions for the transversal component of the velocity and total pressure at the interface between 
the prominence and coronal plasmas, we derive a dispersion relation whose imaginary part describes the evolution of 
the instability. Results are obtained in the limit of weak dissipation. 

Results. Using simple analytical methods, we show that dissipative instabilities appear for flow speeds that are lower 
than the Kelvin-Helmholtz instability threshold. While viscosity tends to destabilise the plasma, the effect of partial 
ionisation (through the Cowling resistivity) will act towards stabilising the interface. For ionisation degrees closer to 
a neutral gas the interface will be unstable for larger values of equilibrium flow. The same principle is assumed when 
studying the appearance of instability at the interface between prominences and dark plumes. The unstable mode 
appearing in this case has a very small growth rate and dissipative instability cannot explain the appearance of flows 
in plumes. 

Conclusions. The present study improves our understanding of the complexity of dynamical processes at the interface of 
solar prominences and solar corona, and the role partial ionisation can have on the stability of the plasma. Our results 
clearly show that the problem of partial ionisation introduces new aspects of plasma stability with consequences on the 
evolution of solar prominences. 
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1. Introduction 

Solar prominences are among the most enigmatic structures 
in the solar atmosphere whose study is made difficult by 
their complex evolution and the multitude of important ef¬ 
fects appearing in them. Prominences are believed to be of 
chromospheric origin and some of them show a long-term 
stability. When formed, prominences maintain their high 
density and low temperature despite being surrounded by 
the million degree solar corona. Their stability and ther¬ 
mal shielding is provided by the magnetic field. Their im¬ 
portance resides in the recognition that almost 80% of the 
observed coronal mass ejections (CMEs) - believed to drive 
the space weather- have a cold chromospheric core believed 
to originate from a prominence, which is why the study of 
the generation and evolution of prominences is necessary. 

The difficulty in studying prominences arises from the 
complex dynamics occurring in these magnetic features, 
but also because of their intrinsic structure and proper¬ 
ties. Early observations showed that prominences are made 
up of fine structures that are composed of many horizontal, 
thin dark threads (filaments) (e.g. de Jager f959, Kuperus 
and Tanclberg-Hanssen 1967) with average width of 200 km 
and lengths of 3,500-28,000 km (see e.g. the review by Lin 
2010 and references therein). These threads are tracers of 
the magnetic field. 

Recent H a and UV/EUV observations showed that so¬ 
lar prominences are also very dynamic with observed bulk 


flows in the range of 2-35 km s -1 (e.g. Berger et al. 2008). 
In active region prominences, flow velocities seem to be 
higher than in quiescent prominences, even reaching 200 
km s -1 , and some of these high-speed flows are probably 
related to the prominence formation itself. The range of 
observed velocities of filament flows is between 5 and 20 
km s -1 . A particular feature in these observations is the 
presence of counter-streaming flows, i.e. oppositely directed 
flows (Zirker et al., 1998; Lin et ah, 2003). Because of the 
physical conditions of the filament plasma, all these flows 
seem to be field-aligned. For a detailed review of the ob¬ 
served flows in solar prominences see Labrosse et al. (2010) 
and Mackay et al. (2010). 

Significant advancement in the study of prominences 
was made when high-resolution observations of waves, os¬ 
cillations, and flows became available. Scientists were able 
to connect theoretical models with observations through 
seismological techniques in order to derive quantities and 
processes (structure of the magnetic field, transport mech¬ 
anisms acting in prominences, internal structure, etc.) that 
cannot be measured directly or indirectly (for a detailed 
discussion of seismological techniques and results see the 
review by Arregui et al. 2012). There is also some evidence 
that velocity oscillations are more easily detected at the 
edges of prominences or where the material seems fainter, 
while they are sometimes harder to detect at the promi¬ 
nence main body (Tsubaki and Takeuchi, 1986; Tsubaki et 
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al., 1988; Suematsu et al., 1990; Thompson and Schmieder, 
1991; Terradas et al., 2002). 

One of the fundamental properties of solar prominences 
is that because of their relatively low temperature, the 
plasma is not fully ionised and its description therefore 
needs special attention, especially when compared to the 
fully ionised coronal plasma that surrounds prominences. 
The ionisation degree of prominences is not well known, 
but there is plentiful evidence that this cannot be ne¬ 
glected when one studies the dynamics and stability of these 
structures (Patsourakos and Vial 2002, Gilbert et al. 2007, 
Labrosse et a. 2010, Zaqarashvili et al. 2011, Khomenko 
and Collados 2012, etc.). 

The problem of prominence stability is paramount for 
other effects such as CME eruption due to the connection 
between these two solar atmospheric structures. In a recent 
series of papers Ryuota et al. (2010), Berger et al. (2010), 
Terradas et al. (2012) highlighted a number physical pro¬ 
cesses taking place in solar prominences that can be con¬ 
nected to instabilities, such as Rayleigh-Taylor instabilities 
(RTI) and Kelvin-Helmholtz instabilities (KHI) under the 
effect of plasma flows. The effect of partial ionisation on 
the stability of prominences was investigated earlier by e.g. 
Diaz et al. (2012), who analysed the appearance of RTI in 
partially ionised prominence plasma. These authors found 
that the linear growth rate is lowered by both the compress¬ 
ibility of the gas and ion-neutral collisions, even though the 
appearance threshold of this instability is not altered. They 
also found that the ion-neutral collisions have a strong im¬ 
pact on the RTI growth rate, which can be decreased by 
an order of magnitude compared to the case corresponding 
to the collisionless limit. They conclude that their results 
could explain the existence of prominence fine structure 
with lifetimes of the order of 30 minutes, a duration that 
classical theories cannot explain. 

In the same year, Soler et al. (2012) investigated the 
KHI of compressional and partially ionised prominence 
plasma. They considered the stability of an interface sep¬ 
arating two partially ionised plasmas in the presence of a 
shear flow. In the incompressible limit the KHI was present 
for any value of the flow, regardless the degree of ionisa¬ 
tion. When extended to a compressible limit, the instability 
threshold was very much sensitive to the collision frequency 
and density contrast between the two layers of their model. 
In particular the density contrast is an important param¬ 
eter in their model. In classical theories the flow speed at 
which the KHI is set is always super-Alfvenic; however, the 
results of these authors show that for a high density con¬ 
trast the threshold can be even sub-Alfvenic thanks to the 
ion-neutral coupling. 

In addition to these instabilities there is another, rather 
unexpected instability that can arise at the interface be¬ 
tween two media called dissipative instability and it is 
strongly connected to the phenomenon of negative energy 
waves. This instability always occurs for flows lower than 
the KHI value. Under normal conditions the interface be¬ 
tween two media allows the propagation of two modes trav¬ 
elling in opposite directions. For flow speeds larger than a 
critical value, the propagation direction of the two waves 
becomes identical, and the wave whose phase speed is 
smaller becomes a negative energy wave (Ryutova 1988). 
The dissipative mechanisms acting in the two regions can 
amplify this negative energy mode leading to dissipative 
instability, and the growth rate of this instability is propor¬ 


tional to the combination of dissipative coefficients. Under 
solar conditions the problem of negative energy waves 
has been studied by many authors (e.g. Ruderman and 
Goossens 1995, Ruderman et al. 1996, Joarder et al. 1997, 
Terra-Homem et al. 2003, etc). In the present study we con¬ 
sider this problem, but the two regions separated by the 
interface are the viscous corona and the partially ionised 
prominence plasma. 

The concept of negative energy waves is based on the 
energy equation 


where E is the linear part of the energy and D is the dissi¬ 
pative function. The two functions appearing in the above 
relation depend on the choice of the frame of reference. If 
we choose a coordinate system where D > 0, then the vari¬ 
ation of the energy with time is negative, meaning that the 
energy of the system decays as a result of dissipation. In 
this case E > 0 for positive energy waves, and dissipation 
leads to the damping of the wave, i.e. to a decay in its am¬ 
plitude. However, if E < 0 the wave is called a negative 
energy wave and dissipation leads to an amplification of 
the wave amplitude resulting in an instability. 

The paper is organised as follows. In Sect. 2 we intro¬ 
duce the basic assumptions regarding the nature of partial 
ionisation and discuss the dissipative mechanisms applied in 
our study. Here we also introduce the mathematical frame¬ 
work we employ when studying the stability of the plasma. 
In section 3 we derive the dispersion relation of incompress¬ 
ible waves propagating at the interface between the viscous 
coronal plasma and partially ionised prominence and es¬ 
tablish the stability thresholds of this model. Finally, our 
results are summarised in section 4. 

2. Governing equations and basic assumptions 

2.1. Equilibrium 

We assume two semi-infinite layers of collisional and incom¬ 
pressible plasma separated by an interface modelling the 
interface between the solar prominence and solar corona. 

The interface between the corona (labelled with index 
“ 1 ”) and the solar prominence (labelled with index “ 2 ”) is 
situated at z = 0 in a two-dimensional (x — z) cartesian 
reference system. The homogeneous magnetic held in both 
regions is along the a:-axis with Hoi 7 ^ -® 02 - The unper¬ 
turbed state is characterised by a magnetohydrodynamics 
(MHD) tangential discontinuity at z = 0, and all equilib¬ 
rium quantities are constant at both sides of the disconti¬ 
nuity. We assume that there is an equilibrium flow in the 
positive x direction in the prominence (for the z > 0 re¬ 
gion), while in the corona (corresponding to 2 < 0) the 
equilibrium is static. The above equilibrium describes the 
interface between a prominence and the surrounding quiet 
corona. Although these two solar regions have been neigh¬ 
bours for a very long time (in the case of quiescent promi¬ 
nences their stability is shown to be of the order of several 
months) they present a very different set of physical pa¬ 
rameters describing them. Solar prominences are cool and 
dense plasma material, thought to be of chromospheric ori¬ 
gin that are surrounded by the very hot and very tenuous 
solar corona. Accordingly, it is customary to consider that 
the density of the prominence is two orders of magnitude 
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larger than the density of the corona and the temperature 
two orders of magnitude lower than the coronal tempera¬ 
ture. Gravity is neglected and so the RTI is not present in 
our problem. 


2.2. Basic assumptions 


The determination of transport mechanisms acting in solar 
plasmas is a very difficult task. After all, the dominant dis¬ 
sipative mechanism depends not only on the location where 
the dynamics occurs, but also on the nature of the phys¬ 
ical mechanisms that needs describing. Under prominence 
conditions Ballai (2003) and Carbonell et al. (2004) showed 
that none of the classical dissipative processes (assuming a 
fully ionised plasma) are able to describe realistic damping 
of observed waves in prominences, except thermal conduc¬ 
tion. Recent studies (e.g. Khodachenko et al. 2004, Arber 
et al. 2007, Forteza et al. 2007, 2008) also showed that 
the dominant transport mechanism in solar prominences 
is probably due to the partially ionised character of the 
plasma. Soler et al. (2009b) found that resonant absorption 
is dominant over ion-neutral effects in the damping of the 
kink mode in prominence threads. In the present study the 
appearance of resonant absorption is prevented by assum¬ 
ing a sharp transition between the prominence and corona. 

In partially ionised plasmas the classical Coulomb re¬ 
sistivity is several orders of magnitude smaller than the 
Cowling resistivity, and the viscosity of the plasma is pro¬ 
vided by the friction between various particles making up 
the plasma (neutrals, ions, protons). The second consider¬ 
ation also implies that the dynamics in solar prominences 
has to be described in a multi-fluid plasma. However, if the 
resistivity of the plasma is dominant (as is assumed here) 
the plasma is described within the framework of single-fluid 
MHD. In the present paper we will assume that these re¬ 
strictive conditions are satisfied, i.e. we are going to use a 
single-fluid description. 

In partially ionised prominence plasmas the Coulomb 
resistivity is many orders of magnitude smaller than the 
Cowling anisotropic resistivity (see e.g. Cowling 1957, 
Khodachenko et al. 2004). Indeed, their difference is given 
by 


Vc = V 


^ 0 [ 


(i) 


where pc is the Cowling resistivity, p is the classical 
Coulomb resistivity, po is the permeability of free space, 
is fraction of neutrals, and the frictional coefficient a n 
in the case of a plasma assumed to be composed entirely of 
H is given by 


Ot-n 


2£n(l-£„) — 


m v 



The degree of ionisation can be characterised by the ioni¬ 
sation fraction (defined as the mean atomic weight, i.e. the 
average mass per particle in units of proton mass) given as 


P = 


2 — 


( 3 ) 


According to this definition, a fully ionised gas corresponds 
to p = 0.5, while a neutral gas is described by p = 1. 

Our aim here is to study the appearance and evolu¬ 
tion of instabilities at the interface of two media, therefore 
we neglect the effects of particle ionisation and recombi¬ 
nation in the solar prominence. Here we assume a strong 
thermal coupling between the species, which means elec¬ 
trons, ions, and neutrals have the same temperature (i.e. 
T e = Ti = T n = T). Therefore, the three-component gas 
can be considered as a single fluid. The concept of a three- 
component gas mixture will introduce new types of trans¬ 
port mechanisms whose importance in the context of solar 
prominences was discussed in detail in the pioneering work 
of Forteza et al. (2007). Since we are going to limit our¬ 
selves to linear dissipation, we will neglect effects connected 
to the inertia of different particles, but also the transversal 
drift of charged particles due to the Hall term and consider 
that thermodynamic quantities (pressure, temperature) are 
relatively smooth functions of the spatial coordinates, i.e. 
the relative densities of neutrals and ions are constants. 
Therefore, when describing the dynamics in solar promi¬ 
nences we will restrict our model to transport mechanisms 
that arise in the induction equation. 

Temperatures in the solar corona can reach millions of 
degrees K, therefore the plasma can be considered to be 
completely ionised. In this important solar region the prod¬ 
uct LOdTi 1 (where ui C i is the ion cyclotron frequency and 
Ti is the ion mean collisional time), therefore ions can gyrate 
many times around magnetic field lines between collisions. 
Under typical coronal conditions this product is of the or¬ 
der of 10 5 . Provided the characteristic scales are larger than 
the mean free path of ions, viscosity in the solar corona 
is mainly due to ions and the viscosity gyrating around 
the magnetic field is given by the Braginskii stress ten¬ 
sor (Braginskii 1965) whose linearised expression takes the 
form of a sum of five terms each with different physical 
meaning (see e.g. Erdelyi & Goossens 2004; Ruderman et 
al. 2000, Mocanu et al. 2008). Under coronal conditions 
the first term, called parallel or compressional viscosity, is 
dominant (by several orders of magnitude) and controls the 
variation along magnetic field lines of the velocity compo¬ 
nent parallel to field lines. The parallel viscosity is due to 
the collision-induced random-walk diffusion of particles and 
is given by 


7 7o = 


PoToksTi 


Too 


where m p is the proton (ion) mass, p is density, ks is the 
Boltzmann constant, and £j „ w5x 10 -15 cm~ 2 is the ion- 
neutral collisional cross-section. The number densities of 
electrons and ions are assumed to be approximately equal. 
The quantity £„ plays an important role in our discussion 
as it contains information about the ionisation degree of 
the plasma. By definition this quantity reflects the number 
of neutrals in the gas mixture, i.e. 

( 2 ) 

p rii + n n 


where po and To are the density and temperature of the 
medium. In practice it is more convenient to work with 
the kinematic coefficient of viscosity defined as u = po/po- 
As determined by Ruderman et al. (1996), a property of 
the highly anisotropic viscosity is that it allows a jump 
in the velocity across a magnetic surface, since a strong 
magnetic field causes ions to rotate around the magnetic 
field lines, thus preventing the diffusion of particles across 
the field lines. This also implies that there is no momentum 
transport across the magnetic surfaces, and different layers 
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of plasma can move with respect to each other without 
friction. 

Observations show that quiescent prominences are made 
of chromospheric material and they live in a relatively sta¬ 
ble position for a long time. High-resolution observations 
very often show that the edges of prominences are not still; 
small- and large-scale features appear and disappear rather 
frequently. We will try to explain these modifications in 
the interface between the two media by instabilities that 
develop owing to the amplification of waves propagating 
along the interface. 

2.3. Governing linearised equations 

We will consider a two-layer system, where an interface sep¬ 
arates the solar prominence and solar corona. The equations 
describing the dynamics of the plasma are the incompress¬ 
ible dissipative and linear MHD equations. In both regions 
the equations 


as suggested in an earlier investigation by Roberts (1981). 
According to his results, in the absence of any equilibrium 
flow, incompressible Alfvenic waves can propagate with a 
phase speed that lies between the Alfven speeds in the two 
regions, which is given by 

^ I ( PlVAl +P2V 2 A2 V /2 | ( V Al+ dv A2 \ 1/2 /, ^ 

k x l, Pl +P2 J { 1 + d J ’ 1 j 

where d = P 2 IPi is the density contrast parameter, k x is the 
parallel component of the wavevector to the interface, and 
vai = Boi/^/popi and VA 2 = B 02 /^/p 0 p 2 are the Alfven 
speeds in the two regions. This dispersion relation describes 
the propagation of the two waves along the interface in 
opposite directions. 

For a stable interface at z = 0 we have to impose the 
linearised kinetic boundary condition and the condition of 
the continuity of normal component of stresses. If v, = 
( v X i , 0, v Z i), then these conditions read 


V ■ v = 0, V ■ b = 0, 


(4) 


are valid. In the solar prominence we assume a field-aligned 
equilibrium flow (vo). As a result, the momentum equation 
becomes 


P 2 - 


<9 v 2 

~dt 


■ v 0 - 


dw 2 

dx 


-VU 2 + 


Bq 2 3b 2 
Ho dx 


(5) 


In the solar corona the equilibrium is static, but the mo¬ 
mentum equation is supplemented by the viscous force, i.e. 



-VPi + 


Bq2 dbi 
Ho dx 


+ V. 


( 6 ) 


In the solar prominence the dominant dissipative effect is 
the Cowling resistivity, therefore the induction equation be¬ 
comes 


3b 2 

~dt 


db 2 

v 0 -r— = B, 
dx 


02 


<9 v 2 

dx 


n. 


(7) 


""“a? = at + ,, 0 a? (11) 

and 

A 

Pi + pivb ■ V(b ■ v) = Pi - P\v-^- = P 2 - (12) 

dz 

We note here that at the tangential discontinuity used in 
the present paper there is no mass transfer between the two 
media, meaning that the state of the plasma in each region 
is not disturbed by the presence of the other medium. The 
system of equations 0 © together with the boundary 
conditions at the interface will form the starting equations 
for our discussion on dissipative instability generated at the 
interface between the two media. 

3. Dispersion relation of surface waves at the 
discontinuity 


In the corona, the flow and dissipative effects do not modify 
the induction equation, so we can write 


3b 1 _ <9vi 

~m ~ Boi lh' 


( 8 ) 


In the above equations vi and b, (i = 1,2) are the per¬ 
turbations of velocity and magnetic field, Pi are the total 
pressures (the sum of kinetic and magnetic pressures), and 
the dissipative terms in Eqs. © and 0 are given by (see 
e.g. Ruderman et al. 1996) 


V = v 


b(b-V)-iv 


b-V(b vi) , 


1Z = ??V 2 b 2 + 


(vc - V) 

|Bo | 2 


V x {[(V x b 2 ) x B 0 ] x B 0 }, 


(9) 


where b is the unit vector in the direction of the magnetic 
field, i.e. b = B 0 /Po- 

Because of the orientation of the equilibrium magnetic 
field the interface between the corona and solar prominence 
can be considered a tangential discontinuity. We write the 
equation of the perturbed interface as 2 = ({x,t). We as¬ 
sume that at |x| — > 00 and |.z| — > 00 all perturbations van¬ 
ish. At the interface surface waves will be able to propagate, 


Since we are going to deal with an eigenvalue problem we 
will perform a normal mode analysis and take all perturba¬ 
tions proportional to exp [i(k x x — uit)] 1 where w is a complex 
frequency that can be written asw = Re(u}) + ilm(u>). This 
particular form of perturbations reduces the boundary con¬ 
ditions to 


Vzi = ~iw(, v z2 = 


(13) 


where fl = ui — k x v 0 is the Doppler-shifted frequency of 
waves in the solar prominence. 

When computing the components of the resistive terms 
given by Eq. © together with the solenoidal condition 
© we can obtain that all dissipative terms containing the 
classical Coulomb resistivity cancel, therefore the dissipa¬ 
tion in the partially ionised prominence is described by the 
Cowling resistivity alone. 

We introduce the viscous and resistive Reynolds num¬ 
bers as 


R v 


vai 

k x v' 


R r 


VA2 
kxVc " 


(14) 


Under coronal and prominence conditions both Reynolds 
numbers are very large and therefore waves will propagate 
with little damping over a period, meaning that in our 
subsequent calculations we will consider that \Re{ur)\ 
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\Im(co)\. The very large Reynolds numbers also allow us to 
consider dissipative terms much smaller than other terms 
belonging to ideal MHD, meaning that in our calculations 
all terms containing v 2 or r \q are neglected. The interac¬ 
tion between flows and waves propagating at the interface 
between the two media is ensured by dissipation that could 
play the role of energy sink. Later we will see that, contrary 
to our expectations, dissipation does not always act towards 
decreasing the wave amplitude; for specific values of flows 
or ionisation degree, the interplay between flows, dissipa¬ 
tion, and waves could lead to an increase of the waves’ 
amplitude, i.e. unstable behaviour. 

In region 1 the viscous MHD equations can be reduced 
to a system of coupled equations for the normal component 
of the velocity vector v z i and total pressure Pi of the form 


dv zl _ _ ik x oj _ 

dz pi (Dai + 2ivk 2 uj) 11 


(15) 


where Da 2 = fl 2 — k^,v\ 2 . Eliminating the total pressure 
from the above two expressions we obtain an equation for 
v z 2 valid in the solar prominence of the form 

= o. (22) 

Using Eqs. m and m we can write that the total 
pressure at the prominence evaluated at the interface be¬ 
haves as 

P 2 = 


P 2 {Da2 + irjcftkpsic 

+ irjck 2 


-^(D^n + iky^rjc). (23) 


The expressions of the two total pressures in the two re¬ 
gions are inserted in Eq. ED. which leads to the dispersion 
relation of the form 



iuco d 2 ^ 

vZid^J 


iu dP± 
Pi'Dai dz ’ 


(16) 


D(oj) = D r + iDi = 0, 
where 


(24) 


where T>ai = w 2 — k 2 v\\- Taking into account that R v 1, 
we can eliminate the total pressure from these two equa¬ 
tions to arrive at a single relation for v z i , i.e. 


D r — Dai + dDA 2 , 


Dt = vk 2 uj + 


dkxVcv 2 A2 

O 


(25) 


d 2 Vzi 7 2 ( i 3 ivklu \ 

dz 2 Dai ) 


v zl = 0. 


(17) 


It is easy to see that v z i will vanish as z —> — oo. In order to 
use the boundary conditions ED and ED we will also need 
to find the value of the total pressure. In order to calculate 
its expression we write the equation for the z-component of 
the velocity ED as 


d 2 v z i 
dz 2 


- a 2 v z i = 0, 


where 


(18) 


a = k x 1 - 


3 ivk 2 u\ 1 
Dai 



3ivk 2 uj\ 

2D A i J ' 


Equation ED allows us to explicitly find the expression of 
the ^-component of the velocity in the solar corona. With 
the help of Eqs. ED and ED we can find that the total 
pressure in region 1 estimated at the interface between the 
two regions can be written as 


In deriving the dispersion relation (1241) we took into account 
that a perturbation method is used meaning that terms 
proportional to v 2 and rfo are neglected. 


3.1. Instability conditions 

Since we assumed that the damping of waves propagating 
along the interface is weak, we can write the frequency of 
waves asw = Re(cu)+ilm(u]) with |Pe(w)| » \Im(u>)\. This 
assumption is in line with our previous statement regarding 
the high Reynolds numbers of solar plasmas and the work¬ 
ing supposition that terms containing squares and products 
of dissipative coefficients can be neglected. According to the 
dependence of perturbations on the variable t assumed ear¬ 
lier, Im(oj) > 0 corresponds to an overstability, i.e. to a 
situation where the amplitude of waves propagating along 
the interface grows as exp (luR). Following the method de¬ 
veloped by Cairns (1979) we write the dispersion relation 
as 


D r (u , k x ) = —iuk^uj — 


idkjr]cVA2 

n 


(26) 


Pi 


PiDai A _ ivujk x \ 
k x \ 2Dai ) 


(19) 


For the prominence we now have an equilibrium flow 
in the positive x direction. Considering again the equations 
that relate the ^-component of the velocity vector and total 
pressure we obtain the systems 



iilcktv 2 A 2 ) , _ Al dP 2 

f 2 J ~ 2 p 2 dz 


and 



irickjv 2 A2 \ dvz2 
f l J dz 



( 20 ) 


( 21 ) 


The solution of the equation D r = 0 will result in the real 
part of the frequency u>. In ideal MHD the interface between 
the two regions is always stable; however, the introduction 
of dissipation may lead to instability. The dispersion rela¬ 
tion for the ideal case can be easily solved and leads to the 
frequency 


Re(uj) ± 


k x vpd 

1 + d 


1 + d . 


d(v 2 KH ~ v o) 


1 1/2 


(27) 


Equation (1271) describes two waves propagating along the 
interface in opposite directions. The quantity vkh is the 
Kelvin-Helmholtz (KH) threshold velocity given by 


v 2 kh = 


1 + d 
d 


(vAi+dv 2 A2 ), 


(28) 


5 

























Ballai et al.: Dissipative instability in partially ionised prominence plasmas 


and it plays a very important role in the discussion of stabil¬ 
ity of waves propagating in a flowing plasma. It is obvious 
from Eq. m that the plasma becomes KH unstable for 
flows that satisfy the condition Vg > v^ H . We estimate the 
value of vkh for the system under investigation. If we con¬ 
sider typical coronal and prominence values for density and 
Alfven speeds (d = 100, vai = 315 km s _1 , VA 2 = 28 km 
s -1 , all taken from Joarder and Roberts 1992) we obtain 
vkh = 423 km s -1 . It is obvious that observations in the 
solar prominences do not show equilibrium flows that are 
larger than vkh ; in reality, these speeds are more likely to 
be of the order of a tenth of vkh or smaller. This means 
that under prominence conditions the plasma at the in¬ 
terface between the solar corona and solar prominences is 
always KH stable. In the absence of any flow, the two so¬ 
lutions of Eq. (1271) describe two modes propagating in the 


opposite direction with equal speeds 


VKH\fd /(1 -I- d) 


In 


the presence of a flow (for our problem the flow is present in 
the prominence while the coronal plasma is at rest), waves 
are drifted by the flow. Since the flow direction points in 
the positive direction, the flow affects the two waves in a 
different way and the symmetry of the two modes is lost. It 
can be easily shown that the difference between the phase 
speeds of the two waves is 2vgd/(l + d). For flow speeds 


larger than vkh/\/ l 1 + d) the direction of the wave prop¬ 
agating in the negative direction is inverted and the back¬ 
ward mode becomes the forward mode. The two modes can 
amplify each other leading to instability. If we plot the two 
frequencies obtained in Eq. ( 1271 ) with respect to an increas¬ 
ing flow speed we obtain that the KHI occurs when the 
oscillation frequencies of the forward and backward propa¬ 
gating surface modes merge for increasing flow velocity. The 
merging point then indicates the threshold of KHI for the 
single interface. In the present analysis we consider flows 
that are less than the KH threshold. 

Since we assume that the damping is weak, the imagi¬ 
nary part of the frequency can be given by 


Im(u )) ~ 


k l ( , dk 2 x ricVA2 \ 

dD r /dRe{u) v n r /'* 


(29) 


where f l r = Re(co) — k x vo- 

Using Eqs. m and 03, it is straightforward to show 
that the imaginary part of the frequencies are 


Im{u>y = 


2(1 + d) 



d{d + 1 )klv 2 A2 ’g c 

2('CoT - dT 2 ) 


(30) 


and 


Im(uj) 


vk l ( v od 
2(i + 0 V r ) 


d(d + l)k 2 v 2 A2 r]c , , 

2O 0 r + dT 2 ) ’ 1 j 


where T = \J d{y 2 KH — Vq). With the help of Eqs. CD)-© 
we can write the Cowling resistivity as 


v 2 A2 m n (2n - 1) f ttto p ^ 1/2 

” C = 2„ 2 (1-„)S,„ ((M ) ■ (32) 

As a result, the two values for the imaginary part of the 
frequency become 


Im[ojy = — 


vk 2 


( vod 
2(l + d) V~ 


+ 1 + 


lm(to) 



Fig. 1. Variation of the imaginary part of the frequency for 
the backward propagating wave with the flow speed and 
the ionisation fraction. The flow changes in the interval 10- 
60 km s _1 and the ionisation fraction varies between 0.5 
(fully ionised plasma) and 1 (neutral gas). The horizontal 
curve is drawn at the Im(ui) = 0 and helps to visualise 
the transition of Im{io)~ from the positive to the negative 
domain. 


did + l)k 2 m n {2iL — l)v\ 2 (nm v \ 

4(z> 0 r-r 2 )(i- M )p 2 E m (lyr) 


and 


vk 2 


2(1 yd) 


v 0 d 


1/2 


-1 - 


d(d + l)k 2 m n (2fi — l)v\ 2 /7r Trip's. 1 ^ 2 

4(u 0 r + r 2 )(l- M )p 2 E m \k^f) 


(33) 


(34) 


We now discuss the sign of these frequencies for a range 
of flow speeds changing in the interval 10-60 km s^ 1 and for 
an ionisation degree varied between the cases correspond¬ 
ing to full ionisation (/x = 0.5) and neutral plasma (/i = 1). 
A simple graphical analysis clearly shows that for the spec¬ 
trum of flows considered here and for any ionisation degree 
the imaginary part of the surface waves that propagates 
in the positive direction (i.e. in the direction of the flow) is 
negative leading to a classical physical damping. In contrast 
the imaginary part of the wave propagating backward (in 
the negative direction) can become positive for flow speeds 
larger than 48 km s -1 (see Fig. 1). A positive imaginary 
part of the frequency is connected to an instability. A con¬ 
tour plot of the imaginary part of the frequency for the 
backward propagating wave is shown in Fig. 2, where the 
role of the partial ionisation and plasma flows becomes ev¬ 
ident. The region above the 3(w) = 0 curve corresponds 
to the region where the wave is unstable, while in the re¬ 
gion beneath the curve the wave is stable and damped. It is 
clear that the flow will destabilise the interface; for a given 
value of ionisation fraction there is a flow value at which 
the interface becomes unstable (a similar conclusion can be 
drawn from earlier studies by e.g. Ruderman and Goossens 
1995). The variation of the zero-level with respect to the 
ionisation fraction shows that as the plasma becomes more 
dominated by neutrals, the plasma interface becomes more 
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0.5 0.6 0.7 0.8 0.9 1,0 


Fig. 2. Contour plot of the variation of the imaginary 
part of the frequency for backward propagating waves. The 
region below the zero level curve corresponds to a stable 
regime and waves will have a classical damping, while the 
interface described by the quantities in the region above 
the curve is unstable. 


stable, so that for an ionisation degree of 0.93 the interface 
becomes stable and waves will damp owing to dissipation. 
Figure |T] also shows that the presence of neutrals stabilises 
the plasma as the instability sets in for higher values of 
flows (here with a density ratio of 100, k x = 5 x 10~ 6 nr -1 , 
v = 10 10 m 2 s -1 , T = 10 4 K, p 2 = 5 x 10~ n Kg m" 3 ). It 
is instructive to identify the role of each dissipative process 
in the appearance of instability. While the partial ionisa¬ 
tion in the solar prominence has the effect of stabilisation 
of the interface, the viscosity in the solar corona will desta¬ 
bilise the discontinuity and the value of the flow at which 
waves become unstable has little variation with the ioni¬ 
sation fraction and significant dependence can be observed 
for larger values of /i. We note that the unstable behaviour 
of the backward wave is also connected to the very high 
density contrast between the solar prominence and corona. 
For a density contrast of one order of magnitude the un¬ 
stable backward wave becomes stable and the imaginary 
part of the dispersion relation describes classical physical 
damping. 

Finally, we explore the connectivity between the dissipa¬ 
tive instability discussed earlier and negative energy waves. 
As specified in the Introduction, the term of negative energy 
wave refers to the situation when the wave energy decreases 
with the increase of the wave amplitude. The energy of a 
wave with amplitude A averaged over one wavelength can 
be given as 




where D r is the dispersion relation of the wave. In this 
case the energy of the wave is the phase-averaged difference 
between the energy of the system when the wave is present, 
and its energy when the wave is absent. A criterion that can 
be used to determine the nature of waves is the formula 
suggested by Cairns (1979) where a wave is considered to 
have negative energy if the quantity 


C = Re(u) 


dD r 

duj 


< 0 . 


(35) 


The function D r is undetermined up to a multiplicative 
constant whose sign has to be determined from the condi¬ 
tion that in the absence of any flows in the system C > 0. 
Comparing this with Eq. (iSl). it is obvious that the condi¬ 
tion for the appearance of dissipative instability is identical 
with the condition of negative energy wave generation be¬ 
cause the expression 

vRe{a? + dklncv^ReW) 
is always positive. 

Another possibility for exploiting the effect of partial 
ionisation on the stability at a magnetic interface is to 
model the interface between two partially ionised plasmas 
of prominences and dark plumes. Observations by Berger 
et al. (2010) revealed the existence of dark plumes within 
the prominence showing turbulent upflows in prominences 
of the order of 15-30 km s -1 . These upflows are believed to 
generate instabilities. In the Ca II H-line plumes are seen 
dark in contrast to the prominence material, which sug¬ 
gests that the plasma in the plumes is hotter and probably 
less dense than the prominence material. The width of the 
plumes ranges between 0.5 Mm and 6 Mm, and their max¬ 
imum heights are between HMm and 17 Mm. The typical 
plume lifetime is between 400 s and 890 s. 

Considering the same prominence/plume parameters as 
in Soler et al. (2012), we obtain that the interface between 
these two partially ionised media becomes unstable for al¬ 
most all values of flows (below the KH threshold) for an 
ionisation degree of the prominence larger than the ion¬ 
isation degree of the plume, but the growth rate of this 
instability is very low, meaning that the dissipative insta¬ 
bility (at least in this simplified framework) cannot explain 
the generation of upflows in plumes by instability. 

4. Conclusions 

In the present study we explored the stability of a tangen¬ 
tial discontinuity by modelling the interface between the 
viscous and fully ionised coronal plasma and the partially 
ionised solar prominence in which the dominant dissipative 
effect is the Cowling resistivity. The magnetic fluids were 
assumed to be incompressible and the prominence equilib¬ 
rium was considered to be dynamical, with a homogeneous 
flow parallel to the interface. Assuming a weak damping 
(confirmed by the very large Reynolds numbers) we ob¬ 
tained the dispersion relation of Alfvenic waves propagat¬ 
ing along the interface. The presence of dissipative effects 
on both sides of the interface renders the dispersion rela¬ 
tion to be complex with the imaginary part of this quantity 
describing the decay or the growth of waves. Our results 
show that while the forward propagating wave is always 
stable, with the amplitude of the wave decaying because of 
dissipation, for the backward propagating wave there is a 
threshold of the flow (below the KHI threshold) for which 
the wave becomes unstable. A careful analysis proves that 
the partial ionisation has a stabilising effect on the interface 
for any value of the ionisation fraction and the unstable be¬ 
haviour can be connected to the viscous nature of the coro¬ 
nal plasma. We also showed that the partial ionisation has 
little effect on the threshold where waves become unstable. 
For a plasma where neutrals are abundant, the instability 
appears for higher values of flows, i.e. neutrals have a sta¬ 
bilising effect. The above results were obtained under the 
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strict restriction of incompressible plasma and a sharp tan¬ 
gential discontinuity between the two plasma layers. The 
same model was used to study the generation of dissipative 
instability at the interface of two partially ionised plasmas 
modelling the prominence and dark plumes. The unstable 
mode obtained in this case shows a very low growth rate, 
meaning that this type of instability (at least in this sim¬ 
plified model) cannot explain the appearance of turbulent 
upflows in plumes that can be attributed to instability. The 
problem of compressibility and a smooth transition between 
the dynamical solar prominence and static corona in the 
presence of ion-neutral friction will be addressed in a forth¬ 
coming study. 
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